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We reexamine the process 7 — > e + +e~ in a background magnetic field comparable to B c = m\je. 
This process is known to be non-perturbative in the magnetic-field strength. However, it can be 
shown that the moments of the above pair production width are proportional to the derivatives 
of photon polarization function at the zero energy, which is perturbative in B. Hence, the pair- 
production width can be easily obtained from the latter by the inverse Mellin transform. The 
implications of our approach are discussed. 

PACS numbers: 12.20.Ds, 11.55.Fv 

The electroweak phenomena associated with an intensive background magnetic field are rather rich. Under a 
background magnetic field, a physical photon can decay into an e + e~ pair or split into two photons. Such processes are 
relevant to the attenuation of gamma-rays from pulsars [0J^]. The study of pair production process 7 + -B — * e + e~ +B 
was initiated by Toll J3] long time ago. He obtained a rather tedious expression for the absorption coefficient «im , 
where || and _L denote the photon-polarization directions which are respectively parallel and perpendicular to the 
plane spanned by the magnetic field B and the photon momentum q. Writing 
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where to is the photon energy and 9 is the angle between the magnetic-field direction and the direction of photon 
propagation. We note that, due to the quantization of electron and positron orbits in the magnetic field, Ku± should 
contain sawtooth absorption edges. However, for B -C B c = m 2 e /e and tosnvO ^> 2m e , these absorption edges are 
rather closely spaced. Hence it is sensible to define an averaged absorption coefficient which is precisely the Tim 
(with a trivial prefactor) displayed above, fn other words, Toll's result is valid for B <C B c and uis'mO 2m e . It 
is interesting to note that, as pointed out by Toll, the functions Tim can not be calculated order by order in eB. 
It is essential to use the exact Dirac wave functions for electrons and positrons in the magnetic field such that the 
resulting Tim are nonvanishing. The non-analytic behaviors of Tim at eB = can be easily seen from its asymptotic 
expression for A <C 1: 

The expression for Tn _l was simplified considerably in the work by Tsai and Erber The authors computed the 
photon polarization function by the proper-time technique |J and determined the absorption coefficient «ii j_ using 
the optical theorem. They obtained 
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where K 2 /3 is the modified Bessel function. At the first glance, the result of Tsai and Erber appears very different 
from Toll's result. However, by computing the moments of Tim , the former authors were able to show that their 
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result is in fact equivalent to that of TollQ. We observe that Tsai and Erber simply utilized the moments of as 
a mathematical tool to show the equivalence between two sets of absorption coefficients. The physical significance 
of these moments was not studied. In this note, we shall clarify the meaning of these moments and develop a new 
method of computing the absorption coefficients. Since our approach is essentially a systematic expansion in B/B c , 
it will remain valid for a background magnetic field comparable to B c . 

We are motivated by the following contour integral, which resembles to the contour integral encountered in the 
QCD sum rule calculation of e + e~ — > hadrons M: 



div 2 II| U (^ 2 ) 
C 2ni (lo 2 + to%) r 
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where the contour of integration C is shown in Fig. 1. 




FIG. 1. The integration contour for and the analytic structure of IIii j_. In actual calculations, we take the radius of the 
circle to infinity. 

The functions H\\,±_ are defined as 

% ± = $ t± n M vei ± , (7) 

where and are respectively the photon polarization vectors parallel and perpendicular to the plane spanned by 
the photon momentum q and the magnetic field B. We note that the integral I n may be evaluated in two different 



To state it more precisely, Tsai and Erber computed the moments of the averaged function T = 1/2 ■ (Tii + T±), which is 
relevant to the attenuation of unpolarized photons. 
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ways. One computes /„ either by the residue theorem or by a direct integration along the contour C with the 
realization that the contribution from the outer circle vanishes. The equivalence of two integration procedures gives 
rise to the relation: 
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where Mn± are the threshold energies of pair productions HQ] given by 
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Mifsin'0 = 4mj, Mian' = m£ 1 + W1 + 2— , (9) 



with # the angle between the photon momentum and the magnetic field. Since k\\,± = lmll\\ ,±/oj by the optical 
theorem, the above equation relate the real part of vacuum polarization function to the absorption coefficient. 

We observe that the l.h.s. of Eq. (0) can be easily calculated at uj 2 — —tu% — 0, since the threshold behaviors of 
Him are absent at this energy value. With this choice of to 2 , we recast Eq. (m into 

1 / d n \ M,If" /-i , , N 
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with y = Mi 2 j_/cj 2 . One notes that the imaginary part of IX || .^(w 2 ) vanishes for the range < uj 2 < M» ± . This 
property has been verified in the previous works |3|,^| . Therefore one can effectively set the integration range of Eq. 
( |To| ) as from y — to y = oo. Now, it is easily seen that the derivatives of IIim at the zero energy are proportional 
to the Mellin transform of k.|| j_ ■ y^ 1 ^ 2 = ■ uj/M^j_. Once the l.h.s. of Eq. ( |Io| ) is calculated, the absorption 
coefficients kim can be determined by the inverse Mellin transform. 

To calculate Him and their derivatives, we begin with the proper-time representation of vacuum polarization 
function Ii^ u in a background magnetic field ||: 
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To construct Him from n^, we note that only the structures proportional to iVn and N± contribute to ELi j_. Since 
we only concern with the limit lo -C ?n e and B < B c , H\\.± can be expanded in a series^] of lo and B: 



2 In fact, we do not need Eq. (|ll]) to obtain such an expansion. A convenient weak-field expansion technique applicable to the 
current problem has been developed in Ref. M. 
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where the neglected terms are of the order (co 2 B 2 sin 2 9/m 2 e B 2 ) n (B/B c ) 2 . Taking the derivatives of IIii ± , we arrive 
at 
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Combining the above equation and Eq. (|T^), the absorption coefficients can be written in terms of inverse 
Mcllin transform: 
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where a is any real number greater than 1/3; while A = (wsin 6B/V3m e B c ) and A" = A' ■ (1 + ^/l + 2B/B c )/2. At 
this moment, we only concern with the leading magnetic-field effect to the absorption coefficients, hence we may set 
(1 + yl-j- 2B/B C ) — ► 2 and A — > A in n±. Numerically we find no distinctions between our results and the results 
of Tsai and Erber 0. For A' = 10, we have K|| = 7.2 x (am 2 /7uj) while kj_ = 4.4 x (am 2 /7rw). For A' = 100, the 
above absorption coefficients become 42 x (am 2 /-7ru;) and 28 x (am 2 /7rw) respectively. In the current approximation, 
km is always greater than n±. For a high-energy photon, i.e., A' 3> 1, we have Km/k_l = 1.5. This is already reflected 
in the above case with A' = 100. For the low-energy photon, A' <C 1, we find Ki|/rtj_ = 2. The numerical agreement 
between Eq. ( |l6| ) and the result of Ref. ||, as shown in Eq. (^), is not a coincidence. We shall verify shortly that 
both expressions arc equivalent by comparing their infinite sequences of moments. 

As mentioned earlier, in order to establish the equivalence of their result with those of the previous works, the authors 
of Ref. § computed the moments c£ ] = J™ d XX n T( X ) with T( X ) = 1/2 ■ {T\\{x) + T ± ( X )) and X = {4m e B c /3ujB). 
The superscript (1) is used to denote the set of moments computed from the functions Tn± given in Ref. §. The 
superscript (0) will then be used for denoting the moments computed from our results for T\\ t ±_- Without taking the 

average, we obtain the moments of each individual function Ty and T±, which we denote as C ( n\\\) and C^\±) 
respectively. We obtain 

,(!),,, n 9 - 2rq , m/2 (6n + 14, 3n + 8) T(n/2 + 2/3)r(n/2 + 4/3)r(n + 2) 
C n >(\\,±) = 2 (3/tt) 1> + 5/2) ' 

To compare Cn (||, J- ) with the moments pertinent to Eq. (|i"o|), i.e., Dn = J °° dy ■ y"" 1 • (T\\ ,±{y)y~ 1 l 2 ), we employ 
the relation 

where C , 2°' ) _2 (II' -L) is defined in the same way as c|^_ 2 (|| , -L). Our first goal is to show that C^_ 2 ( \\ , _L) = C2„_ 2 (||,-L). 
In fact, this identity can be established by combining the relation 

n(0)ni | s _ 2 2 " +1 ( B \ 2n - 1 r(3n - l)r(2n)r(2n) (6n + l,3n+l) 
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derived from Eqs. (|10|) and (|15|), with Eqs. (|17j), (|1J) and the identity 
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Now that we have shown C^_ 2 (||, J-) = C$_ 2 (\\, JL), we obtain the identity £>4 0) (||,-L) = D&^dl, J-) where D^dl.-L) 
is given by the r.h.s. of Eq. (O) with C 2 °_ 2 (||, -L) replaced by C^^dlt J-). Since both Z?„'s are identical, one can 
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show that the absorption coefficient derived from our approach, Eq. ([l6|), is equivalent to the result of Ref. || given 
by Eqs. (||) and (1), provided that 
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according to the Carleman's theorem jl(|. Indeed, this is true since 
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Therefore, we have proven that the leading-order results of our approach agree with the results of Ref. ||. In addition, 
as one can see from the r.h.s. of Eq. (|l8|), only the even moments defined by Tsai and Erber are relevant to the 
physics of pair-production in a background magnetic field. 

We like to point out the differences between our approach and the approach of Ref. |5| . Tsai and Erber begin with 
IV^v given in Eq. ( pi] ) and evaluate the imaginary part of for uj = q° greater than the pair production threshold. 
They arrive at the asymptotic result, Eq. (^), in the limit B <C B c and wsinf? ^> 2m e . However, their approach 
does not provide an estimate of possible corrections as B and uj deviate from the above limit. Our approach has an 
advantage in that it treats the magnetic-field effects perturbatively for B < B c . In Eq. ([l5]), the n-th derivative of 
the vacuum-polarization function Hu j_ is expanded in powers of B 2 / B 2 . Hence the absorption coefficient, which is 
related to the derivatives of II|| __l by an inverse Mellin transform, can also be written in powers of B 2 /B 2 . In this 
way, we are able to compute the absorption coefficient even for B comparable to B c . As for the low energy regime 
near the pair production threshold, wsin# > 2m e , the quantum effects due to the magnetic field become important. 
Namely, for given uj and B, the momenta of e + and e~ along the magnetic-field direction can only take discrete values, 
and consequently the absorption coefficients K\\ j_ contain resonant peaks. The spacing of these peaks increases as 
wsin6* gradually decreases to the pair-production threshold 2m e . A detailed study of this threshold behavior has been 
initiated by Daugherty and Harding pj{ ]. In Fig. 6 of Ref. pT[ , it is shown that the threshold behavior is non- negligible 
for £ = uo 2 B c /2m 2 B < 10 3 with sin 2 9 = 1. For a general 9, the relevant parameter becomes £ = uj 2 sin 2 9B c /2m 2 B. 
It should be understood that our result as well as the result of Ref. f| are applicable for a large £(£ ) where the 
threshold effect is not significant. 

We have mentioned that our results for Kii j_ are written as power series in (B/B c ) 2 . It is important to compute 
the next-to-leading corrections. Let us begin by computing the next-to-leading magnetic-field corrections to n || and 
its derivatives. To do this we perform a weak-field expansion with respect to the exponent 4>o and the amplitude , 
along with a rotation of the integration contour s — > —is: 
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Hence the next-to-leading (NL) corrections to the derivatives of nil reads: 
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Using Eq. (|To|), and applying the inverse Mellin transform, we arrive at 
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where b can be chosen to be any positive number. Numerically, for (B j B c ) 2 = 0.1 and A' 



10, we have n¥ L 
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! x (am 2 /iruj). We note that £ w 10 4 for the current values of B and A . In this case |ty /«|| | < 0.1%. For the same 
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magnetic-field strength with £ = 10 3 (A' ~ 3), we obtained «|| = 2.0 x (am 2 /7ra;) and «j| L = 7.7 x 10 x (am 2 /iru)). 



The ratio r = | /"c/X ^ /« n | remains to be less than 0.1% in this case. Hence in the energy regime that the quantum 



effects of the magnetic field is not essential, the subleading contribution to the absorption coefficient, is rather 

suppressed. If we extrapolate our analysis down to the energy of pair-production threshold uj sin 8 = 2m e while 
maintaining (B/B c ) 2 = 0.1, i.e., A' = 0.35, we find K|| = 6.5 x 10" 3 x (aml/nuj) and nf L = -2.2 x 10~ 3 x (aml/iru;). 

It is interesting to see that is of the same order of magnitude as the leading contribution. This reflects the 
limitation of our approach and that of Tsai and Erber near the pair-production threshold. 

The next-to-leading correction to n± is calculated in a similar way. The expansion of (f>o proceeds as before while 



cosh^v) v cosh(z) sinh(zw) 2(cosh(zw) — cosh(z)) 
sinh(z) sinh 2 (z) sinh 3 (z) 

= - JL(1 - t- 2 )(3 + v 2 ) + ^(1 - v 2 )(15 - 6v 2 t, 4 ) + • • • . (26) 



Then the next-to-leading corrections to the derivatives of II_l are 
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Applying the inverse Mellin transform, we obtain 
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where c can be chosen as any positive number, and, to isolate the 0(B / B 2 ) corrections, we have made the identifica- 
tions A -> A and + + 2B/B C ) -> 2. Numerically, for (B/B c ) 2 = 0.1 and A' = 10, we have n± — 4.4 x (am 2 /iru>) 
while = 1.3 x 10~ 3 x (aml/nui). Similar to the K\\ case, \k± l /kx\ < 0.1% for A' = 10. For A' = 3 which corre- 
sponds to £ ps 10 3 , we find «_l = 1.2 x (aTO 2 /7Uj) and = — 4 x 10~ 3 x (am 2 /7rw). In this case \k i [ l /k±\ ps 0.3%. 
We observe again that the ratio, \kt[ /«ll|> grows rapidly to 60% at the energy of pair-production threshold with 
k± l being negative. 

^From the above next-to-leading order calculations, it is quite evident that the 0(B 2 / B 2 )- corrections to kii and 
k± are both rather insignificant. However, one should be reminded that there are still 0(B/B c )- corrections to k± 
as shown in Eq. jl^) . Without making the identifications A — » A and (1 + y/T+ 2B/B C ) — > 2, we have k± = 3.8 x 
(am 2 e /Tiuj) for A" = 10 x (l + ^/l + 2B/B c )/2 (i.e., A' = 10) with (B/B c ) 2 = 0.1. We recall that k± = 4.4 x (am 2 /7rw) 
if the above identifications are made. Hence the 0(i?/B c )-correction reduces k± by about 14% at the current u> and B. 
Essentially, the 0(_B/i3 c )-corrections to k± relative to K|| are kinematic in nature. They are due to the differences in 
pair-production threshold resulting from different polarization states of the decaying photons. For a photon with the 
1 1 -polarization, the pair-production threshold corresponds to both e + and e~ being in the ground state. On the other 
hand, for a photon with the _L-polarization, cither e + or e~ must be in the first excited state at the pair-production 
threshold [|. 

Since the 0(i? 2 /-B 2 )-corrections to K\\ t ± can be neglected, Eq. ( |l6| ) are accurate expressions for photon absorption 
coefficients even to a magnetic-field strength comparable to B c . In Fig. 2, we plot K\\ x as functions of u> with 
sin 2 8 = 1 and (B/B c ) 2 =0.1. Here our result for km is identical to that of Tsai and Erber, whereas our result for k± 
contains the kinematic corrections which were not taken into accounts in previous works. As can be seen from Fig. 2, 
both «|| and k± are rather insensitive to lo once lo surpasses 10 MeV. Hence, for a sufficiently large ui, the kinematic 

corrections to n± mainly reside in the factor 1/(1 + y/l + 2B/B C ) rather than in A . It should be noted that the 
applicability of our result is determined by the parameter £ = to 2 sin 2 8B c /2m 2 B |ll|. For the parameter set in Fig. 
2, i.e., (B/B c ) 2 = 0.1 and sin 2 8 = 1, we have lo ps 25m e = 12.5 MeV for £ = 10 3 . Hence, one expects the plots in 
Fig. 2 to be reliable for u > 10 MeV. 
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FIG. 2. The photon absorption coefficients as functions of the photon energy with sin 2 6 = 1 and (B/B c ) 2 = 0.1. 

At this point, one might conclude that our approach can not describe the threshold behavior of the absorption 
coefficient. This is in fact not true. We should stress that it is the weak-field expansion of Eq. ( |l5| ) that spoils the 
threshold behavior of the absorption coefficient, despite such an expansion is useful for computing the absorption 
coefficient at much higher energies. To see this, it is instructive to review how Eq. ( |l5| ) is derived. In this note, 
we obtained Eq. (|l5|) by expanding Eq. ( |ll| ) directly. However, one could also derive Eq. ([l5]) by expanding 
the internal fermion propagators in Him in powers of eB. This expansion has been derived in Eq. (47) of Ref. 
P|. In that equation, one can see that the weak-field expansion mixes contributions coming from different Landau 
levels. Hence the threshold behaviors of Him at any given Landau level are spoiled by the weak-field expansion. In 
fact, the absorption coefficients km j. calculated from Hii ± do not vanish below the pair-production threshold, i.e., 
< lu 2 < Afj?j_. This can be seen by performing the integration in Eq. (|l6|). Such a behavior is again an artifact 
caused by the weak field expansion performed near the pair-production threshold. Indeed, as mentioned before, ftj^x 
are comparable to K\\ ± with an opposite sign in this energy range. The cancellation of these two contributions is 
consistent with the fact that the absorption coefficients should vanish below the pair-production threshold 

To recover the threshold behavior, one may compute n Mt , using electron propagators in the Furry picture [I12l[bJ : 
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i.S n \u> ,q 2 ,q 3 )a/3 j 



(29) 



where q 2 = (q ) 2 — (q 3 ) 2 , n is the quantum number for the Landau level, and S n - iUtq 2 q 3 is a 4 x 4 matrix in the 
spinor space. In this form, all the poles of the propagator appear explicitly, and the threshold behavior of is 
preserved throughout the calculation |14|] . Clearly, a calculation using the Furry-picture propagators compliments 
the weak-field expansion technique we have been discussing so far. The former produces a correct threshold-behavior 
of the absorption coefficient but becomes unpractical at larger energies, since, in such a case, contributions from a 
great number of Landau levels has to be summed over. Nevertheless, in the scenario that B 3> B c , it is convenient to 
use the Furry-picture propagators because the available Landau levels for pair production to occur are significantly 
reduced. 

In conclusion, we have developed a new method for computing the photon absorption coefficient in a strong back- 
ground magnetic field B < B c . Disregarding the next-to-leading magnetic-field corrections, our k\\ is identical to that 
obtained by Tsai and Erber jjj. Although Tsai and Erber derived Ku under the assumption B <C B c = m 2 /e and 
ui sin 6 3> 2m e , we have been able to show that such a result is in fact accurate for B comparable to B c , provided u> 
(£ ) is large enough. For k±, our result differs from that obtained in Ref. Q. In this regard, we have identified certain 
0(B I B c ) corrections to k±_ which are kinematic in nature. We also pointed out that our approach may be extended 
to lower photon energies near the pair-production threshold, so long as we calculate n M „ with electron propagators 
in the Furry picture For a supercritical magnetic field B 3> B c , we argued that it is convenient to use the 

Furry-picture propagators. As a closing, we like to emphasize that a better understanding of the current process is 
crucial for determining the photon attenuation properties in highly magnetized pulsars 
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